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We analyze the weak gravity in the braneworld model proposed by Dvali-Gabadadze-Porrati, 
in which the unperturbed background spacetime is given by five dimensional Minkowski bulk with 
a brane which has the induced Einstein Hilbert term. This model has a critical length scale r c . 
Naively, we expect that the four dimensional general relativity (4D GR) is approximately recovered 
at the scale below r c . However, the simple linear perturbation does not work in this regime. Only 
recently the mechanism to recover 4D GR was clarified under the restriction to spherically symmetric 
configurations, and the leading correction to 4D GR was derived. Here, we develop an alternative 
formulation which can handle more general perturbations. We also generalize the model by adding 
bulk cosmological constant and the brane tension. 
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A braneworld model, whose gravity behaves as four dimensional at short distance scale but it shows higher dimen- 
sional nature at larger distances, was proposed by Dvali, Gabadadze, and Porrati (DGP) [1]. In this model, the brane, 
on which the fields of the standard model are confined, has the induced Einstein- Hilbert term [2,3]. This model has 
various cosmologically interesting features [4-10]. Particularly in the model with five dimensional bulk, an interesting 
cosmological solution was found, in which accelerated expansion of the universe at late epoch is realized without 
introducing the cosmological constant [7]. Based on this model, a novel mechanism that dilutes the cosmological 
£NJ . constant was also proposed [10]. 

Although we have mentioned above that the gravity in this model at short distances is expected to behave as 
four dimensional, it is not so transparent if the model actually mimics four dimensional general relativity (4D GR). 
The linear analysis of this model shows that the tensor structure of the induced metric perturbations takes the five 
dimensional form even at short distance [1]. The situation is analogous to the case of models with massive gravitons. 
In this case the deviation from 4D GR does not vanish even in the massless limit, which is known as the van Dam- 
Vcltoman-Zakharov discontinuity [11,12]. In this context, the possibility that the 4D GR is recovered by non-linear 
effect was suggested by [13]. There have been many discussions about this issue [14]. We have in particular a 
' clear statement that the discontinuity disappears when we introduce the cosmological constant [15,16]. Although the 
analysis with the cosmological constant is quite suggestive, the discontinuity is absent only when the limit is taken 
keeping the length scale determined by the cosmological constant much smaller than the compton wavelength of the 
j^Tj, massive graviton. The length scale determined by the cosmological constant must be longer than the Hubble horizon 
size. Hence, the recovery of 4D GR cannot be proven by introducing negligibly small cosmological constant as far as 
. £^ j the graviton mass is not completely negligibly small. 

Also specialized to the context of five dimensional DGP model, there are various works aiming at answering the 
question whether 4D GR is recovered at short distances or not, and many evidences that indicate the recovery of 
4D GR were reported [17-21]. It was shown that the evolution equation for homogeneous isotropic universe becomes 
identical to that for 4D GR when the Hubble expansion rate is much larger than the inverse of the critical length 
scale, r" 1 [17]. In Ref. [19] it was clearly shown that the linear analysis breaks down at scale shorter than (r^r^) 1 / 3 
since the brane bending becomes non-linear there. Further, a consistent form of a black hole metric induced on the 
brane was presented. Approximate black hole solution including the bulk was constructed under the restriction to 
spherically symmetric configurations [20]. The same paper also gave the leading order correction to 4D GR, which 
is potentially observable by the future development of precision measurement of solar system [22,21]. The results in 
Ref. [20] were extended to the case with the background of expanding universe [21]. 

In this paper, we develop an alternative formalism which can handle general perturbations in weak gravity regime. 
To handle general perturbations, we restrict our consideration to the case that the unperturbed metric on the brane is 
given by Minkowski spacetime. We also make further generalization to the model that also takes into account the bulk 
cosmological constant and the brane tension balanced with it. (Such generalized DGP model was discussed before in 
Refs. [23-26].) We confirm the recovery of 4D GR at short distances and rederive the leading order correction to it. 
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II. SETUP 



The model that we consider is defined by the five dimensional action 

s=^|rf 5 ^(fl (5) +^)+ j d i x ^^^ R w- 3 M +Lmatter y (2 .i) 

where M^,r c and I are constants. R^ and R^ are, respectively, the curvature scalars corresponding to the five 
dimensional metric g M „ and the four dimensional one g$ induced on the brane. Here, we added both the bulk 
cosmological constant and the brane tension terms to the original DGP model. They are tuned to admit the Minkowski 
brane as a vacuum solution. The model is reduced to the original one by setting I — > oo. The unperturbed background 
geometry is given by five dimensional anti-de Sitter spacetime 

ds 2 = g [ ^dx a dx h = dy 2 + 7 ^{y)dx p dx v = dy 2 + e- 2y /% v dx»dx u , (2.2) 
with a brane located at y = 0, where Z 2 -symmetry is imposed. Here rj^ is 4 dimensional Minkowski metric. 



III. SEMI NON-LINEAR PERTURBATIONS 



We follow the method of Ref. [27] introduced for the purpose of analyzing weak gravity in the Randall-Sundrum 
model [28]. We prepare two coordinate systems. In the coordinates {x a }, the gauge is chosen so that the metric 
perturbations h ab can be easily computed in the five dimensional bulk. Namely, we use the Randall-Sundrum gauge, 



h 5a = 0, 



V = 0, 



0. 



(3.1) 



In this paper the fifth direction is the direction of extra-dimension. The greek and latin indices represent four and 
five dimensional coordinates, respectively. The other coordinate system {x a } satisfies the Gaussian normal conditions 



h 5a = 0, 



(3.2) 



and also keeps the location of the brane unperturbed at y = 0. Under the coordinate transformation x a = x a — <; a (x), 
the metric perturbation transforms as 



h ab = h ab (x - £(*)) + ( -g ( ° b ] 5 e + ^9 { Z^ 5 ) 2 ■ 



- {e, a (g%Hx -0 + h cb (x £))) + (a ~ b)} + (g%\x - £) + M* " 0) 



(3.3) 



The argument of the variables is supposed to be x unless otherwise is specified, and " )0 " denotes a differentiation with 
respect to x a . 

The conditions that the {00} component and {0/i} components are zero in both coordinates provide equations for 
the gauge parameters, which are solved up to second order as 



(2) 

where £ 5 (x p ) and £^(x p ) are the values of the gauge parameters evaluated on the brane, and 

Jo 

= f ' dy-f" 
Jo 

We assume the following order counting, 



(3.4) 



(2) 



(2) 

S> ,pS ,<7 S ,/ 



f 5 f b - 

£ 5 < e — — f p h p I 



<e 2 , 



(3.5) 



(3.6) 
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and keep the terms up to 0(e 2 ). Here e 2 is the order of the Newton potential $ = — ^/ioo- Later we will verify the 
consistency of this order counting. Then, the transformation for {\iv\ components reduces to 



h^(x) = h^(x - £(x)) + dhftv, 



(3.7) 



with 



(3.8) 



Hereafter, the Greek indices are lowered or raised by the metric -f^ u . 

The brane location is given by y — 0. Hence in the {x a }-coordinates the brane is bent. For simplicity, we impose 
the harmonic gauge condition, v — for the induced metric on the brane. To the second order in e, this 

condition gives 



□ " 



7& + (°?)£ 



(3.9) 



From this relation, we find that the assumption ^ v — 0(e 2 ) is consistent if the assumed order of £ 5 is correct. 
Substituting (3.9), the gauge transformation Sh^ v evaluated on the brane becomes 



Then, the trace of the induced metric is also evaluated as 

h = + 2D" 1 ( 7 ^VTL£W - (°f) 2 ) • 

Next we consider the junction condition. After a straightforward calculation, we can show 



(d g + 2l~ l ) V = (d y + 2T 1 ) h„ v {x - + + , (at y = 0), 



with 



Using this relation, the junction condition becomes 



Ml 
2r c 



2S„ 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



The equation that determines the brane bending is obtained from the trace of the above equation as 



1 



h 



cn 1 — ttT = 

Ml 2 r c 



e + □■ 



-l " 



(3.15) 



where we have introduced r* := [(l/r c ) + (2/.£)] _1 . We can neglect the last term on the right hand side of Eq. (3.15) 
because it is always higher order compared with the first term. The left hand side is something like the Newton 
potential <&, hence we assume it to be 0(e 2 ). Outside the matter distribution with the total mass m, the left hand 
side can be expressed as w r g /r, where r g := to/47tM|. At large scale, r > (^g^ 2 ) 1 ^ 3 , the first term on the right hand 
side dominates, while at small scale, r < (rgr 2 ) 1 / 3 , the second term dominates. Therefore we have 



C>(e 2 ) = max 



(3.16) 



Thus we find that our assumption as to the order counting for £ 5 is justified. 
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IV. MECHANISM FOR RECOVERING FOUR DIMENSIONAL GENERAL RELATIVITY 



The remaining task is to evaluate (dy + 2^~ 1 ) h^x — £_(x))\y =0 . Here we need to solve the bulk field equations. 
Different from the R-S case, we solve the bulk equations with the Dirichlet boundary condition (3.7). Here, we note 
that the location of the brane is not a straight sheet in the coordinates in the R-S gauge {x a }. 

We can give the general solution for the bulk field equations as a superposition of homogeneous mode solutions 
with purely outgoing boundary condition: 



K„(x) = J H^(py P ^K 2 (p£ey/ e )d A p 



(4.1) 



where K 2 (p£) is the modified Bessel function and is the expansion coefficient. The coefficient TL^ U is to be 
determined so as to satisfy the Dirichlet boundary condition 



Vb=o(*) = / H^py^-^Kzipie^^d^p 



(4.2) 



If we are allowed to approximate the above expression by setting £ a = 0, we have h fll/ (p) := 
(27r) -4 J d 4 x e~ %v f xP h^ix) = 'H ) _ lv (p) 1 and therefore we have 



K 2 (pl) 



K 2 (pt) 



(4.3) 



We think that the errors caused by this naive approximation are not large, although any rigorous proof is not ready 
yet. If the leading errors are simply proportional to h^, we can neglect them since they are higher order in e. Such 
a naive expansion with respect to £ will be justified for small p. But for large p, we will not be allowed to expand 
the combination p£ in the exponent. However, as we will see below, even the leading correction to the gravitational 
potential coming from the contribution of this part is suppressed to be irrelevantly small at small scale r < r c . Hence, 
the errors due to this naive approximation can be crucial only if this approximation significantly underestimate the 
magnitude of (dy + 2£~ 1 )h lll/ \y =0l which is quite unlikely. 

Using (4.3), the junction condition (3.14) is written down explicitly as 



V 



( p 2 2pK 1 (p£) \ as 2 ~~ 



r c K 2 ( P £) 



(at y = 0), 



(4.4) 



where 



and 



V 



, , V K Y (p£) 



p z + 



r c K 2 (pl) 



6h$ = n- 1 W°£ 5 J% U - 2(nf )i%„ 



(4.5) 



(4.6) 



Here the quantitiy with ~ represents the Fourier coefficient of the corresponding variable as before. We show that in 
the square brackets on the right hand side of Eq. (4.4) the first term gives the dominant contribution. We can drop 
the last two terms simply because they are always higher order in e compared with the first term. The second term 
is irrelevant since it can be eliminated by a four-dimensional gauge transformation. As a result, the equation that 
determines the metric induced on the brane is reduced to the one for the linear theory The only difference is in 
the equation that determines the brane bending (3.15). 

The order of magnitude of the first term on the right hand side of Eq. (3.15) is estimated as £ 5 /r c <~ $ 1 / 2 r/r c at 
small scale, and hence it is suppressed by the factor <I>~ 1//2 r/r c compared with the Newton potential $. The leading 
behavior of the induced metric is therefore determined by setting the left hand side of Eq. (4.4) to zero. Thus we 
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conclude that 4D GR is recovered by taking into account the non-linear brane bending for weak gravity at small 
scale r < (r^rg) 1 ^ 3 . If we take the limit r c — > oo, all length scales come into this regime. Hence, we have confirmed 
the absence of van Dam-Veltman-Zakharov discontinuity. However, because of the factor $ -1 / 2 the leading order 
deviation from 4D GR at small scale is less strongly suppressed than the naively expected suppression r/r c , as first 
pointed out in Ref. [20]. 

At large scale, this term becomes more and more important. For r > (r^r g ) 1//3 , we have 

(4.7) 



2„2 ' 



3M 4 V 

This is nothing but the result for the linearized case. In the next section, we discuss the regime where the linear 
theory is valid. After that, in the succeeding section, we discuss the leading order correction to the 4D GR at short 
distance scale assuming static and spherically symmetric configurations. 



V. LINEAR REGIME 



In this section we consider perturbations at large scale r 3> (r^rg) 1 ^ 3 , where the linear theory is valid. Substituting 
Eq. (4.7) into Eq. (4.4), wc obtain 



M'i 



K" ~ T72 V ( - ^ory^f ) , (5.1) 



with 



a= — - t . (5.2) 

r c + 



r c + 3 



2 



First we look at the behavior of the propagator D, which is already discussed in Ref. [26]. When we consider the 
length scale much smaller than £, we have 

, (^»1). (5.3) 

V + r c P 

At the length scale smaller than r c (pr c 3> 1), the propagator V behaves as that for a four dimensional field. On the 
other hand, at the intermediate scale between r c and I (r c » p^ 1 S> £), the propagator behaves as that for a five 
dimensional field. When the length scale is much larger than I, Ki(pl)/ K 2 (p£) goes to pl/2. Thus we have 

Hence again the propagator V behaves as a four dimensional field, but the Newton's constant is not given by 2M^ 2 
but by 2M 4 - 2 /(i +£/2r c ). 

Next we turn to the tensor structure specified by a. For four dimensional massless graviton we have a = 1, while 
a = | for the case of massive graviton. For pi ^> 1, we have 

(5-5) 

We have a — > 1 for I <C r c , while a — > 2/3 for I <C r c . On the other hand, for pi < 1, we have a — ► 1 irrespective of 
the ratio between t and r c . 

The results are summarized in Fig.l. When r > £, the 4D GR is realized by the Randall- Sundrum mechanism. 
The effective Planck mass differs from M4 in this case. On the other hand, when r < r c , the gravity becomes four 
dimensional again, but the tensor structure differs from 4D GR. 
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FIG. 1. Summary of the results of linear analysis, which is valid when r 2> (rj?r g ) 1//3 . The horizontal axis represents length 
scale. The raws labeled T> and a, respectively, explain the property of the propagator T> and the indicator of the tensor structure 
a. 



VI. STATIC SPHERICAL SYMMETRIC NON-RELATIVISTIC STAR 



We consider a static spherical symmetric non-relativistic star. Assuming that the energy momentum tensor is 
dominated by the {00}-component Too = p, we neglect the effect of pressure when we solve the metric perturbation. 
First we solve non-linear equation for £ 5 (3.15). Under the present assumptions, Eq. (3.15) is simplified as 



i 



d r 



3r k 



e- 2 K£) 2 



This equation can be immediately integrated once, and we obtain 



3r 



where 



1 



r g( r ) 



= 0, 



Outside the star, we have r g (r) = r g = m/AwM^. Solving the above equation with respect to £ 5 r , we have 




(6.1) 



(6.2) 



(6.3) 



(6.4) 



Here we have chosen the signature in front of the square root imposing the condition that £ 5 r does not become large 
at r — > oo. The other branch with "+" sign is outside the scope of the present formalism since we have assumed the 

Minkowski brane background. At small scale, this expression reduces to £ 5 w - 



dr \/r g (r)/2r. Outside the matter 



distribution, we simply have £ 5 = —y / 2rr g . Hence, the correction to the Newton potential is given by 

<5$ p 



(6.5) 



which recovers the result obtained in Ref. [20]. 



VII. CONCLUSION 



We developed a formalism to calculate the metric perturbations induced by the matter localized on the brane in 
the generalized DGP model, in which the bulk cosmological constant and the brane tension terms are added. Here 
we clarified the mechanism for disappearing the van Dam-Veltman-Zakharov discontinuity in this model. In our 
approach, the crucial point was to take into account a part of second order perturbations of the brane bending. This 
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method was largely motivated by the recent works by [29,30]. Our scheme almost completely controls the order of 
magnitude of the neglected higher order correction terms. In this sense, we think that this work gives an alternative 
sufficiently satisfactory proof of the absence of the van Dam-Veltman-Zakharov discontinuity in the DGP model. 
Under the restriction to the static and spherically symmetric source, we confirmed that our formulation correctly 
reproduces the leading order correction to 4D GR at short distances obtained in [20] . 
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